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Abstract 

A new model of gamma bursts (GB) generated by MACHO's in the form of non- 
compact dark matter halo objects (NO) is proposed. The model explain well the 
Q> ' spherical symmetry, logN — logS curve of GB absolutely identically to the well-known 

■ giant dark matter halo (GDMH) model. The main difficulty of GDMH model - the 
^ i lack of neutron stars in our galactic halo is eliminated because of the large number of 
pi ^ NO (~ 10^'^) which is quite close to a total number of bursts during the life time of 

■ the Universe. The proposed mechanism of the burst is the instability and explosion 
, of the baryonic core (BC) of NO. The structure of BC and possible mechanism of BC 

I \ explosion are discussed. 

The intensive infrared emission from BC of neutralino stars (NeS) at the wavelength 
^1 (3— 8)/i is predicted. The possibility to discover NeS by observation of infrared emission 

CD ' from the recently distinguished point-like gamma sources is indicated. 

(N : 

O 1 Introduction 

i> : 

^ , The origin of gamma bursts is one of the mostly intriguing problem in astrophysics during 

'"^1 the last 20 years. Analysis of the observational data shows that today exist two main models 

^ \ of the origin of gamma bursts: cosmological and giant halo model [1]. 

■ Giant dark matter halo (GDMH) model was proposed in [2]. It considers relic neutron 

■ stars as the source of gamma bursts. These stars should have the same density distribution 
• • , as distribution of cold dark matter in the Galactic halo, which has a size Rh ~ 200kpc 

• ^ \ [2,3]. This proposition allowed to explain the main statistical properties of gamma bursts: 

^ \ their tightly spherical symmetry and a considerable condensation to the center, known as 

■ logN — logS curve [1]. Furthermore GDMH model fairly well describes weak asymmetries 
in the distribution of gamma-ray bursts, observed during first years of measurement by the 
BATSE device of COMPTON Observatory [4,5]. The characteristic energy released in 
gamma-rays in one burst as follows from this model is 10^^ — 10^^ erg. 

The weak point of this model is the supposition that the source of gamma bursts is 
relic neutron stars (NS). The full number of these stars is not large enough N < 10^ — 10^, 
what means that every NS must repeat gamma burst many times: 10^ — 10'' , since the 
total number of bursts should be ~ 10^^. The same difficulty has analogous model, which 
consider high velocity NS as a source of gamma bursts [6,7]. These stars have got high 
velocities during supernova explosion and are moving out of the Galaxy. The number of 
such NS is about 10^ — 10''. It is difficult to find the mechanism of repeating gamma bursts 
from NS and the source of energy needed for repeating 10^ times bursts. It should be noted 
also that according to observations bursts never repeat (we do not speak here about special 
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small group of the so-called "repeaters"). These difficulties for the NS models as the source 
of gamma bursts were emphasized in numerous papers [2-9]. 

Here we will propose a new mechanism of gamma burst connected with the instability 
and a possible explosion of the baryonic core of noncompact dark matter halo objects (NO). 
Microlensing established, that the large part of dark matter in Galaxy halo consist of MA- 
CRO'S - unseen objects with masses (0.05 — 0.8)Mq [10,11]. These objects are naturally 
considered to be White Dwarfs, Brown Dwarfs, Jupiters or some other stars or planets. To 
the contrary in [12,13] was supposed the existence of quite a new type of massive structures 
- gravitationally compressed objects from nonbaryonic cold dark matter (CDM). These ob- 
jects are the results of a small scale hierarchical structuring developing in CDM [14,15]. They 
are noncompact, having some spectrum of characteristic dimensions and masses. After re- 
combination a baryonic component falls down into the potential well formed by nonbaryonic 
matter. In this way the compact baryonic core is created. The nonbaryonic component 
forms a spread spherical halo, where the dominant mass of the object is concentrated. The 
main hypothesis of [12]-[15] is that these noncompact nonbaryonic objects (NO) determine 
significant part of microlensing events observed in [11]. 

NO were proposed in [12,13] to avoid the difficulties arising from microlensing observa- 
tions. Here we consider them as the source of gamma bursts. Because of supposition that 
the main part of dark matter consist of NO this mechanism gives the same statistical results 
as GDMH model. It explains well the spherical symmetry, logN — logS curve and small 
asymmetry absolutely identically to [2-5]. The new very essential point is connected with 
Nno - the number of NO in Galaxy halo. According to [12-15] N^o ~ 10^^ what is quite 
close to the total number of bursts in our Galaxy during lifetime of the Universe. So no 
bursts repetition is needed for this model. 



2 Dark matter in the baryonic core. 

Non dissipative cold dark matter creates stationary spherically-symmetrical structures com- 
pressed by their own gravitational forces. For these objects the theory estabhshed a funda- 
mental scaling low for density distribution [16]: 

p oc r~", a fti 1.8 (1) 

which holds quite well for large scale structures in the Universe. 

Noncompact dark matter objects (NO) arc the small scale structures with the masses 
Mx ~ (0.01 — 1)Mq and characteristic dimensions Rx ~ 10^^ — 10^^ cm [13-15]. For small 
scale dark matter structures the distribution (1) holds on up to the dimensions Tc an order 
magnitude less than [17]: 

rc ~ O.lRx- (2) 
The density distribution of CDM particles in NO could be presented in a form [17]: 

< r < Tc 

Vc <r < Rx (3) 



P 
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This distribution does not take into account the existence of baryonic core. The funda- 
mental difference between nonbaryonic dark matter and baryonic is determined by dissipa- 
tive processes. After recombination the baryonic matter under the action of radiative coohng 
loses its thermal energy and falls down into the potential wells formed by nonbaryonic mat- 
ter. It leads to the creation of baryonic body (baryonic core) with mass ~ O.OSMj. at 
the center of the dark matter object. 

It is essential, that the process of the baryonic matter condensation change significantly 
the density and a full amount of the nonbaryonic matter trapped inside the baryonic core 
(BC). Let us consider this process in details. One can neglect the collisions between non- 
baryonic particles, which are oscillating in the potential well created by their gravitational 
field. It is easy to estimate, that the typical oscillation time (~ 10 years) is essentially less 
than the time of baryonic core formation. It means that the core formation is a slow process. 
Because of this reason we can describe the changing of distribution function of dark matter 
particles using the adiabatic approximation. Initial adiabatic invariant /j is defined by the 
formulae: 



Here E = e/rUx, e-is the energy of a particle moving in potential por^, V'oo- is the depth 
of the potential well, m- is the particle angular momentum and Tmini^max are reflection 
points, determined as a roots of the expression under the integral. Since we are interested 
in changing of distribution function near BC, in (4) was chosen (in accordance with (3)) 
V' = -000 + 27rG'po/3r2. 

After the formation of BC the adiabatic invariant // of a particle near the center takes 
a form 



here GMb/r -is a potential of a created baryonic body. In expression (5) we took into account 
only potential of baryonic body and neglected the potential of dark matter. As we will see 
below (12) the mass of CDM particles in BC and its neighborhood is essentially less than 
the mass of baryonic body. 

In view of conservation of the adiabatic invariant, it follows that for any value of /, we 
have for the distribution function /(/): 





(4) 




(5) 






It is natural to choose the initial distribution function /o as a maxwellian one 






(7) 
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After integration of (6), (7) , one can find 



Substituting (8) into (6) we obtain: 



Let us find the asymptotic of solution (9) in the hmit r — > 0. In this case E — oo, and 
taking into account, that < 2GMf,r, we find: 

It is easy to sec that the law for density distribution (10) does not depend on the shape 
of distribution function (7), but it depends on its behavior at — The law (10) has a 
singularity at r — 0. The behavior of density is defined by this singularity of potential, but 
in reality the singular law (10) is cut on the scale of baryonic core ri,. Borrowing in mind (in 
accordance with (3)) that 

^Trnom^r^ ^ O.IM^, (11) 
and choosing T pa ■^qo one can find from (10), (11): 

Pin) = 33 X nom^ — — 

So we see, that the adiabatic capture of dark matter particles lead to the increase of the 
dark matter density inside the baryonic core in 10^ — 10^ times. Because of this the full 
amount of the mass of a dark matter trapped inside the BC Mj;b could reach a significant 
part: 



It means for M^, = 0.5Mq, that M^^b ~ 10^^ - 10^^ g. 



10-6 _ 10-7 (12) 



3 The structure of baryonic core with CDM compo- 
nent. 

Let us consider now the structure of baryonic core. The stationary state of this core is 
described by equations: 

VP = -Vi/j 

V (kVT) = -aT^S + Qn + Qth (13) 
AV' = AnG {pb + Pn) 
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Here P = Pb + Pn is the total pressure of baryons and nonbaryons, T(r) is the temperature 
of baryonic particles, ■0- potential of gravitational field, K-heat conductivity, pb and p„ are 
densities of baryons and nonbaryons. Qn and Qth are the sources of the heating of BC 
by DM particles. These sources arc significant for the DM particles of a special type and 
would be considered in the next section. is the temperature of a core surface S, a - 
Stephan-Boltzman constant. 

The distribution of CDM particles inside the baryonic core is determined by their colli- 
sions with baryons,it means that in stationary conditions they have Boltzman distribution: 



Pn = Pno exp 



rp-k 



(14) 



where ip{r) is gravitational potential and T* is the effective temperature of CDM particles. 



Jo 



Rt 



Ef{E)dE 



'^+f^ ^E-'ilj- m?/2r 
W)dE 



> dr 



(15) 



r^+^ \lE-il)-rn?/2r^ J 
Here p„o is CDM density in the center of the core: 




Equations (13) take into account that the crossection of CDM particles with baryons is small 
enough. It means that these particles oscillate many times inside the core between collisions 
with baryons. Using (15) we obtain CDM term in equations (13) 



dPn 
dr 



T* dp, 



dr 



Equations (13) have two different solutions describing two different stationary states of 
the baryonic core containing nonbaryonic component. The first one is usual, determined by 
baryons in the whole range of r: < r < r^. It has value of baryons density in the center 
p^l^^Q — PbO and distribution Pbir) gradually falling down to the edge of the body p — > 1, 
r — > rb- For the planet of Jupiter type pbo ~ 30 g/cm^, but the temperature inside Jupiter is 
essentially less than in the center of BC. In our case the mass of baryonic body Mb ~ O.OSMr 
could be higher, than Jupiter mass, but taking into account temperature, one can expect, 
that central baryonic density by the order of magnitude is 



PbO ~ 10 



9 



cm 



(16) 



If the density of nonbaryonic component Pn(r) is everywhere less, than Pb{r) 



Pn{r) « Pb{r) 



(17) 
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the existence of CDM particles, trapped inside BC will lead to the small corrections to the 
usual baryonic stationary state. 

But according to the Boltzman distribution (14) 

Pn = Pno exp |- (^^^ I (18) 
Here r„o gives the characteristic scale of the heated region 



/ ST* 

\ 27:0 pborux ^^^^ 

In (18), (19) we supposed a constant baryonic density pf, — p^o ■ Density of CDM particles 
in the centrum of the core p„o is connected with full mass M^b — by relation: 

^"0 ~ 3/2 3 y^^) 



TT 



nO 



Condition (17) is fulfilled only if 



n 



' nO 



< Pw 



Taking into account (19) we see, that for low temperatures T* or for the high mass of CDM 
particles if: 

MnPlJo y gy* ) > 1 (21) 

condition (17) is not fulfilled. In this case p„o > Pbo and CDM particles could significantly 
affect the density distribution being trapped by their own gravitational field in the central 
part of BC. General stationary solution of equation (13) is shown in the Fig.l. One can 
see from the figure that two different stationary modes for potential exist. One is the 
mode determined mostly by baryons, this mode we have discussed earlier. Another is the 
CDM particles selftrapped mode. In the vicinity of the centrum r — > the last mode is 
asymptotically described by: 

IT* / r \ 

^^-—ln{-^) + ... (22) 
P 



TTlx \Rc 



27rma;Gr 

Here Rc determines the dimension of the region of selftrapped CDM bulk 

G Mrfn^ , , 

= i-^ 23 

Mc is the full mass of selftrapped neutralino. The main term in (22) coincides with well- 
known Boltzman distribution of noninteracting gravitationally selftrapped matter [18]. 

The existence of a second stationary state means that if in the process of the slow cooling 
of the core conditions (21) could be reached then the transformation to the second stationary 
state will begin. This reconstruction will lead to the loss of the hydrodynamic stability of 
the system resulting in the release of a significant amount of energy E ~ 10'^"'' — 10"^"^ erg 
(see (35), (36)). This energy will be emitted into a wide spectrum of electromagnrtic waves 
, including gamma rays. X-rays, optic emission and radio waves. 
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4 Explosion of baryonic core in neutralino star 



The particles of which CDM is composed are not known now, though there exist some 
hypothetical candidates: neutralino, heavy neutrino, axions, strings. Neutralino and heavy 
neutrino like the Major ana particles can annihilate in mutual collisions. Noncompact objects 
which consist of such particles are partially disappearing because of annihilation processes. 
These specific objects we call "neutralino stars" (NeS) [12] and for nonbaryonic particles 
from which they consist of we will use as a general name "neutralino" [13]. 

Annihilation of neutralino leads to effective radiation of gamma photons and heating of 
baryonic core. Under definite conditions the core can explode. It should be emphasized, that 
the energy released in all these processes was input into CDM during the period of its freezing 
out. The possibility to detect CDM particles through their annihilation was discussed in a 
number of papers (see [19,20]). In these papers always was supposed a smooth distribution 
of CDM in the Galaxy. The existence of NeS leads to a new situation , connected with 
compression of CDM in NeS, which results in strong amplification of neutralino annihilation 
processes and gives significant constrain on possible type of neutralino particles [13,17,21]. 
Prom the other side it opens the possibility of existence of quite observable fiuxes of gamma 
radiation by NcS. In particular, in [17] a NeS model of diffusive gamma emission and recently 
discovered nonidentified point like gamma sources [22] is developed. Here we will discuss a 
model of gamma bursts, which is connected with the possibility of blowing up of the baryonic 
core of NeS. 

As was shown in section 2, the amount of the trapped in BC neutralino particles is high 
enough M^^ ~ (10^^ — 10^^)Mx. Because of this reason the energy released in neutralino 
annihilation could effectively heat the baryonic core. On the other hand, the neutralino 
density p„ could be small in comparison with the density of baryons pi,. 

The heating term Qn in equation (13) could be written in a form 



where Pn{i^)- is the neutralino density (18)-(20), and r„ is the lifetime of neutralino in baryonic 
core. Taking into account that neutralino annihilation is a rare event in comparison with 
neutralino collisions with baryons, we obtain the following expression for 



Here < av > is a characteristic crossection of neutralino annihilation multiplied on typi- 
cal velocity. Note that t„ does not depend on r. So, the heating efficiency Qnis directly 
proportional to Pn{f)- 

The distribution of Pn{f) is given by equation (18). It follows from (24), that if conditions 
(21) are fulfilled the heating is going mostly in the central region of the core. 

Let us estimate now the main parameters of the core. From (19) it follows that: 



Qn — 




(24) 



2.4 X lO'^ {r^/lQ^K^f^ ilQGev/m;)^''^ (p&o/lO^ 



-1/2 



(26) 



cm 



cm 
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One can see that the dimension of heated region could be not too small in comparison with 
the radius of baryonic core: 

R^^ / 3Mb y^'^2xl0i° cm (27) 



Here and below we take: 

Mb = 0.05M^ - 5 X 10^^ g, = 0.5Mq = 10^^ g, pb =< Pb >= 1 — , 

cm 

Density of neutrahnos in the center of BC, as follows from (20) is: 

Pno- 1.3x10 [-^g)[-^) [y^) [To^^) ^ ^2^) 
and their number density: 

M„, \ flO'^K^V^' r \V2 / V^' 



no = 7.6 X 10^^ 



lO^Vy \ T* J \10GevJ \10g/cm 



We see that conditions (17), is usually fulfilled. 

Characteristic annihilation time for neutralinos in the core could be of the same order as 
the lifetime of the Universe: 

. 3.2 X 10- fi^l C^V" 1'-^] (^^V^ C^V" . (30) 



<av>o J \ T* J \ Mn J \lQg/crrf) \ m. 

Here we took into account that according to [13], the dominant annihilation process of 
neutralinos is p-wave with crossection: 

The annihilation time of neutralino should be greater than life-time of the Universe to- It 
depends very strongly (according to (30)) on the parameters of the NeS. If it is not so, the 
mass of neutralino M„ will decrease due to annihilation until t„ ~ to- It means that energy 
release due to annihilation which heats central region of the core in present time is (24): 

Qno = 6.8 X 10""-^ (32) 
cm s 

Consequently the full energy dissipated by neutralino in BC is: 

Q = Qnor^o^'/' = 1-4 X 10^9 ^ (33) 

s 

This energy being transported to the surface of the core is emitted due to usual blackbody 
radiation. 

<j = g,=,^5=1.4xlO»(^)'(j^)^f^ (34, 
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Relation (33) determines the surface temperature as a function of parameters. For the 
stationary state formulas (25)-(33) describe some average situation, the exact distribution of 
temperature would be established in accordance with equation (13). To fulfill the stationary 
conditions a convective heat transport near the boundary region has to be selforganized. 

Let us estimate the hydrodynamic stability (13) comparing the full pressure of baryonic 
gas P with the gravitational energy E. Taking constant we have: 

£=i^ (35) 

were is the full mass of the baryonic core. Total particle energy 

P^Nf 

where T is the average temperature and N - number of baryonic particles. Here we took into 
account that the temperature is high enough to keep the object in gaseous state. Taking 
f 0.6T* = 6 X 10^^°, one can obtain 

P fti E = 2.5 X 10^^ ( 5 JiO^^g ) ^^^"^ 

We see that the particle energy is quite close to the gravitational one, it means that the 
main integral relation for equation (13) is well fulfilled. 

So the structure of baryonic core of NeS because of neutralino heating is more alike to the 
star than to the planet. In the stationary conditions it has the temperature in the central 
part of the core an order of 10^K°. But the surface temperature is not very high. 

We emphasize that the time needed to reach stationary temperature distribution in BC 
is compatible with the age of Universe. Because of this reason nonstationary processes could 
be significant. As the heating is going in the central part the temperature in this part 
T* at some conditions could be much higher than the stationary one. Nonstationarity of 
T* depends on the thermal conductivity and grows especially strongly with the full mass 
of neutralinos M„ trapped in the core. The strong growth of T* follows from the relation 
(24) for neutralino dissipation rate and if M„ > 10^^ g it could increase on two orders of 
magnitude or even higher. This process can lead to the effective heating of the small central 
part of the core . The temperature in the center of BC could then increase in nonstationary 
process up to 

T^^IO^X". (37) 

So we see that the baryonic core of NeS because of the heating by neutralinos is quite 
peculiar astrophysical object which could be constructed either like a planet or like a star 
depending on the efficiency of heating. The core can lose its hydrodynamical stability and 
explode, releasing the energy of the order of (35) 



~ (0.1 - 0.01) P ~ 2 X (10^^ - 10^^) erg 



This energy is enough for creation of gamma burst {E^ ~ 10 — 10 erg), just what is 
needed for GDMH model of gamma bursts. We will consider now three main destabilization 
factors which can lead to explosion of BC. 
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1. Overheating of the central part of BC by neutrahno. 

A very strong heating of BC is connected with the high value of bulk of trapped neu- 
tralinos M„ > 10^^ g. The heating would be more concentrated in the centrum of BC for 
the high values of neutralino mass m^. As the surface temperature of BC is of the order 
500° — 1000° (what follows from the conservation of thermal flux) it is only weakly ionized 
what leads to the strong depression of thermal conductivity. Convective transport would 
be the main process of the heat transport in this region. Overheating in the central region 
could lead to the loss of hydrodynamical stability and explosion of BC. 

2. Thermonuclear heating. 

Neutralino annihilation inside the star leads to generation of highly energetic protons and 
gamma photons. These energetic particles produce ions of D2 and some other light elements 
by interaction with the bulk protons of the core. Because of this process the number of D2 
and other particles is constantly growing in time in the central part of BC. Following (24), 
one can estimate the number density of deuterium produced for the lifetime of the Universe: 

ND2=lD2-^^to (38) 
rrixC 



ND2-2X10 ^n2[j^^j [j^) [j^^^ ) (,^7To^ 

7£)2 is the transformation coefficient for Gev-energy protons into D2 ions due to collisions with 
the bulk protons. The number density (37) is big enough for initiation D2 He^ He^ 
thermonuclear reaction. It means that if the temperature could reach high values (36) the 
effective nuclear process will begin. This process determine the heating term Qth in equation 
(13) Thermonuclear process leads to a very effective heating Qth (13) of the central part of 
BC and could become explosive , depending on initial state. It is necessary to mention that 
in (37) the total mass of trapped neutralino should be take into account, because even if 
according to (30) netralino annihilate, the process of production of D2 will take place. 
3. Selftrapping of neutralino. 

This process has been already considered in the previous section. As was mentioned 
above the dimension could be small, it is diminishing with the growth of the mass of 
neutralino particles ttIx- On the other hand to reach the conditions (21), which determine 
the possibility of formation of selftrapped neutralino bulk (SNB) easier for low temperatures. 
So, it could be that SNB is formed during nonstationary process, when temperature in BC 
reach its minimum. 

The fundamental feature of SNB is that the gravitational selftcompression is growing 
up with the growth of temperature T*: its dimension became smaller (23) and neutralino 
density becomes higher (22). Such process leads to an explosive heating of the central region 
of BC. As follows from (22), (23), (24), (31) the energy, dissipated in SNB in 1 cm ^ per second 
is: 

= 3.2 X 10- f {—)' f^)' ( 
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And the full dissipated power is: 



Were r^j„ is a minimal dimensions for the bulk. We see, that the energy dissipation in SNB 
is explosive, lasting only a few seconds. 

It is evident, that to reach in reality such a stationary state is impossible , but non- 
stationary transition to this state is also a fast process, which can explode the BC. Let us 
determine the characteristic time scale r„ of this process. Neutrahnos crossection for elastic 
collisions is usually lower than for annihilation. So, they can establish Boltzman distribu- 
tion only through collisions with baryons. Characteristic crossection for neutralino-proton 
collision according to [23] is: 

<.„,„>«io--f « ) ^ 

yig/cm J s 

Supposing the density of protons (16) pb ~ 10 g/cm^ we obtain 

Tn ~ lO^S 

This time is much smaller than the characteristic time of establishing of thermal equihbrium 
by the heat transport. So we conclude, that it is possible to have an overheating and blowing 
up of BC due to selftrapping, heating and fast compression of neutralino bulk in the centrum. 

5 Conclusions 

1. We have shown here that MACHO's in the form of noncompact CDM halo objects 
could serve as the source of gamma bursts. This model could be proved if in experiments 

the significant local asymmetry of gamma bursts distribution, connected with the existence 
of a Giant Halo of the Andromeda (M31) would be established by observations. To do that 
a device with the sensitivity an order of magnitude higher than BATSE is needed. These 
experiments were already proposed earlier in connection with GDMH model [5]. 

2. The explosion of barionic core of NO proposed here as a source of gamma bursts will 
lead also to the emission of X-rays, optic infrared and radio waves. So a wide spectrum of 
electromagnetic emission follows the gamma burst. 

3. Let us discuss now the interesting new opportunities for observations which follow from 
our model. The intensive infrared emission from BC of neutralino star (NeS) is predicted in 
present work according to formulae (34). The wave length of this emission A ~ (3 — 8)//, its 
full power Q ~ 10^^ — 10^^ erg/s. Emission should have blackbody spectrum with effective 
temperatures Tg ~ (500 — 1000)/^°. The intensity of this emission at the Earth from the 
nearest sources should be of the order S ~ (10""'^° — 10^^^) erg/s cm^. 

We emphasize, that the theory predicts gamma emission from the same NeS also. In [17] 
are distinguished nonidentified point-hke gamma sources from EGRET data [22] . It is shown 
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that these sources could be considered as the emission generated by NeS. A more detailed 
analysis [15], demonstrated the existence of two types of the NeS gamma sources. Galactic 
type has the emission intensity an order of magnitude higher, than halo type, considered 
in [17], what is in full agreement with the intensity of point-like gamma sources observed 
by EGRET [22]. In [15] demonstrated also an especially good agreement of a new EGRET 
observations with the predictions of a space distribution of sources made in [17], basing on 
NeS model of gamma sources. So of a great interest is to observe the infrared emission from 
the same sources. If this observations would be successful it would mean that we see the 
neutralino stars: neutralino halo of such a star radiates in gamma and its baryonic core in 
infrared rays. 

We conclude, that infrared observations of nonidentified point-like gamma sources [17] 
are of fundamental significance as they could lead to the discovery of NeS. 

The authors are grateful to V.L.Ginzburg for permanent interest and fruitful discussion, 
and to I.Axford, V.A.Dogiel, A.Lukyanov and V.Sirota for valuable remarks. 

We are grateful for hospitality to the Max Planck Institute of Aeronomy, where a part 
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Fig.l 

Two solutions of equation (13). First one correspond to neutralino selftrapping regime 
(curve 1) at the center of barionic core, and the second one is regular solution defined 
basically by barions density (curve 2). 
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